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ABSTRACT 

We  study  the  stability  of  viscometric  flow  using  the  type  of  short  memory 
introduced  by  Akbay,  Becker,  Krozer  and  Sponagel  [1-7].  The  instability  found 
by  these  researchers  is  recognized  as  a  change  of  type  leading  to  non¬ 
evolutionary  character  of  the  governing  equations. 

We  also  address  the  question  of  justification  for  the  short  memory 
assumption  and  find  that  it  cannot  be  justified  for  some  of  the  more  popular 
rheological  models. 
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SIGNIFICANCE  AND  EXPLANATION 


— In  flows  of  viscoelastic  liquids,  instabilities  frequently  arise  at  high 
shear  rates.  The  explanation  for  these  instabilities  is  still  a  largely 
unsolved  problem.  This  paper  discusses  boom  consequences  of  a  recent  theory 
due  to  E.  Becker  and  his  coworkers  at  Darmstadt.  This  theory  predicts 
instabilities  if  the  first  normal  stress  differences  is  sufficiently  high.  We 
show  that  this  instability  is  associated  with  a  change  of  type  in  the 
governing  equation. 

The  theory  is  based  on  a  certain  approximation.  The  justification  for 
this  approximation  is  also  discussed  and  is  found  to  be  problematic.  _ . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sussMry  lies  with  MPC,  and  not  with  the  authors  of  this  report. 


REMARKS  ON  THE  STABILITY  OF  VISCOMETRIC  FLOW 


M.  Ahrens*,  D.  D.  Joseph**,  M.  Renardy*  and  Y.  Renardy^ 

"Complex  problems  have  simple,  easy-to-understand  wrong  answers." [9] 
1.  INTRODUCTION 

In  viscometric  flow  of  polymeric  liquids,  instabilities  are 
frequently  observed  at  high  shear  rates.  In  many  cases  these 
instabilities  occur  even  th  ough  the  Reynolds  number  iS  low 
and  have  distinctly  different  characteristics  from  the  instabilities 
that  arise  from  inertial  effects.  Although  various  mathematical 
theories  have  been  developed,  the  proper  explanation  of  these 
instabilities  is  still  an  open  problem  (for  review,  see  Petrie 
and  Denn  [18]).  Among  the  mechanisms  proposed  are:  Changes  of 
type  in  the  governing  equations  which  may  lead  to  fracture  or 
phase  changes  in  the  material,  development  of  shocks,  free 
surface  instabilities,  slip  on  the  wall,  or  bifurcations  of  the 
classical  type.  In  general  context,  equations  describing  linear 
perturbations  of  viscometric  flow  were  derived  by  Pipkin  and  Owen 
[19]  and  Dunwoody  and  Joseph  [11] .  These  equations  provide  a 
general  framework;  however,  in  order  to  arrive  at  definite 
conclusions,  more  specific  assumptions  are  necessary.  Such 
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assumptions  can  be  of  two  kinds:  First,  a  particular  constitutive 
model  many  be  assumed,  and  second,  one  may  seek  approximate 
equations  valid  in  certain  limiting  situations.  Either  approach 
seems  to  have  its  problems.  In  the  former  case,  the  analysis  of 
one  constitutive  model  tells  you  nothing  about  others,  and  it 
is  not  clear  at  this  time,  which  models  are  adequate  for  which 
materials — in  particular  at  high  shear  rates.  In  addition, 
the  determination  of  the  stability  boundaries  for  simple  shear 
flow — either  analytically  or  numerically — does  not  appear  to  be 
easy  even  for  the  simplest  constitutive  models  [13] .  On  the 
other  hand,  "approximate"  constitutive  laws  must  be  used  with 
great  caution.  There  may  be  problems  of  consistency  (e.g.  does 
instability  occur  in  a  range  of  parameters  consistent  with  the 
assumptions  needed  to  justify  the  approximation?)  or  stability. 

That  is,  the  stability  properties  of  approximate  equations  need 
not  be  the  same  as  those  of  the  "exact"  equations,  as  is  well 
known  for  "fluids  of  grade  n",  if  n  >  1  [10],  [14],  [21], 

An  approximate  theory  of  perturbed  viscometric  flow  based  on 
a  hypothesis  of  short  memory  (small  relaxation  time)  has  been 
developed  by  U.  Akbay,  E.  Becker,  S.  Krozer  and  S.  Sponagel  in 
Darmstadt.  In  a  series  of  papers,  they  study  various  flow  situations 
such  as  plane  Couette  and  Poiseuille  flow,  flow  between  eccentric 
rotating  cylinders  and  boundary  layer  flow  [1-7] ,  [20] .  The 
purpose  of  this  paper  is  to  work  out  some  consequences  of  this 
theory  and  to  give  a  critical  assessment  of  it.*  The  short 
memory  hypothesis  can  be  formulated  in  terms  of  a  perturbed 
extra  stress  in  the  form 

•The  theory  of  Akbay  et  al.  is  explained  in  Chapters  2  and  3. 
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(1.1)  S(t)  *  j  K ( K, s) D( t-s)  ds 

where  K  is  assumed  to  decay  rapidly  to  zero  fast  and  to  have 
appreciable  values  only  for  s  <<  k"*1.  E.g.  we  could  take 

(1.2)  K(tc,s)  -  A(K)e"s/X 

where  kA  <<  1.  We  have  denoted  by  D  the  symmetric  part  of  the 

perturbed  velocity  gradient,  and  K  is  a  fourth  order  tensor 

so  that  (KD) .  ■  K .  _  D«  .  k  is  the  shear  rate  of  the  basic  flow. 

“  i  j  x  j  xm 

The  approximation  applied  to  (1.1),  assuming  (1.2),  is  that  only 
terms  linear  in  kA  are  kept,  while  higher  powers  are  neglected. 

The  short  memory  hypothesis  is  applied  only  to  the  perturbation 
equations,  not  to  the  "basic"  flow  that  is  being  perturbed.  The 
viscometric  functions,  i.e.  the  viscosity  n(<)  and  the  normal 
stress  differences  N^ftc)  and  N2(k)  are  assumed  to  be  given 
independently  of  any  short  memory  hypothesis.  As  we  shall 
see  later,  this  leads  to  serious  consistency  problems  at  least 
for  some  of  the  constitutive  laws  presently  used  in  rheology. 

The  theory  of  Akbay  et  al.  would  work  best  for  a  constitutive 
law  that  has  a  different  relaxation  time  for  perturbations  of 
simple  shear  flow  than  it  has  for  simple  shear  flow  itself. 

It  is  not  clear  to  us,  how  such  constitutive  laws  should  be 
constructed . 

Akbay  et  al.  make  the  additional  approximation  that  the 
Reynolds  number  is  neglibible.  With  these  assumptions  they 
find  that  the  following  condition  is  necessary  for  the  instability 
of  plane  Couette  flow  to  two  dimensional  disturbances. 


<l-»  [h 

u  K 

A  similar  criterion  is  obtained  for  plane  Poiseuille  flow.  They 
also  give  numerical  evidence  suggesting  that  (1.3) is  sufficient 
for  instability  to  occur.  More  recently  they  have  also  considered 
three  dimensional  disturbances  [20] . 

In  this  paper,  we  reconsider  the  stability  of  plane  Couette 
flow  and  we  show  in  Chapter  5  that,  in  the  limit  of  zero  Reynolds 
number,  (1.3)  is  indeed  a  necessary  and  sufficient  condition  for 
instability.  The  instability  is  associated  with  a  change  of  type, 
a  notion  which  we  explain  in  detail  in  Chapter  4.  The  character 
of  this  instability  is  quite  different  from  those  familiar  from 
Newtonian  fluid  mechanics.  At  the  critical  "Weissenberg  number" 
not  only  one  mode  becomes  unstable,  but  a  whole  family  of  modes. 

At  any  supercritical  Weissenberg  number,  ther  is  an  ®  number  of 
unstable  inodes.  These  unstable  inodes  have  arbitrarily  short 
wavelengths,  and  their  growth  rates  actually  tend  to  »  as  the 
wavelength  tends  to  zero.  As  a  consequence,  this  kind  of 
instability  must  lead  to  a  much  more  severe  change  in  the  flow  field 
than  the  development  of  secondary  flows.  "Fracture"  may  be  one 
possible  consequence  of  such  an  instability. 

For  finite  Reynolds  numbers,  the  criterion  (1.3)  is  still 
valid  for  a  change  of  type.  We  do  not  know  whether  particular 
modes  may  become  unstable  at  lower  Weissenberg  numbers,  but 
numerical  results  indicate  that  this  is  not  so  at  low  Reynolds 
numbers.  In  fact,  for  disturbances  of  a  fixed  streamwise 
wavelength  the  critical  Weissenberg  number  seems  to  increase  with 
the  Reynolds  number. 
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In  a  final  chapter/  (Chapter  6) ,  we  attempt  to  give  a 
critical  assessment  of  the  theory  of  Akbay  et  al.  First/  we  find 
that  there  is  a  problem  with  the  short  memory  assumption.  It 
turns  out  that  if  this  assumption  is  applied  to  the  basic  flow, 
then  the  normal  stresses  are  small  compared  to  the  shear  stress, 
and  (1.3)  cannot  be  satisfied.  We  give  a  general  argument  for 
this,  and  also  illustrate  it  for  a  particular  model.  It  thus 
appears  that  (1.3)  and  the  short  memory  hypothesis  are  inconsistent, 
unless  the  fluid  somehow  has  a  different  memory  for  perturbations 
of  shear  than  it  has  for  shear  flow  itself.  We  do  not  know  if 
and  how  this  could  be  made  precise.  Another  question  that  might 
be  asked  is  whether  the  result  of  Akbay  et  al.  would  hold  without 
the  short  memory  hypothesis.  The  upper  and  lower  converted 


Maxwell  model  satisfy  (1.3)  at  high  enough  shear  rates,  in  fact  the 
quotient 


becomes  arbitrarily  large.  However,  it  was  recently  shown  [16],  [22] 
that  a  change  of  type  does  not  occur  in  those  models.  Thus  any 
existing  instability  would  not  be  of  the  kind  considered  by  Akbay 
et  al.  Whether  there  are  instabilities  at  all,  is  a  difficult 
problem.  Numerical  attempts  to  find  them  have  so  far  led  to 
negative  results  [13] . 


2.  PERTURBATION  OF  THE  STRESS  IN  NEARLY  VISCOMETRIC  FLOW 


In  this  section  we  derive  general  equations  for  the  perturbed 
part  of  the  extra  stress,  with  no  restriction  on  the  memory  of 
the  fluid.  The  essential  features  follow  the  method  of  Dunwoody 
and  Joseph  [11] .  We  refer  the  reader  to  that  paper  for  a  detailed 
derivation. 

Consider  a  viscoelastic  fluid  in  which  the  extra  stress  is 
given  by  a  functional  of  the  deformation  gradients 

(2.1)  3  =  3T  -  T  +  PI  *  ^|"g( s)~l 

-L-  -*  s=o 


Where  T  is  the  Cauchy  stress,  P  is  the  pressure  and  G  is  the 

T 

relative  Cauchy  strain:  G(s)  =  Vx  (s)VX(s)  -  I_,  where  X(s)  = 
X(x,t-s)  is  the  position  at  time  t-s  of  the  particle  presently 
(i.e.  at  time  t)  at  x.  We  represent  the  basic  flow  (as  yet 
unspecified)  by  a  "o"  subscript,  and  write  the  perturbation  as 
follows: 


(a) 

X(s)  = 

X(S) 

+  £(s) 

_o 

(b) 

U  (s)  = 

+  u(s) 

— o 

(c) 

p  =  Po 

+  p 

(d) 

3  =  s 

—  =o 

+  S 

Here  the  soleroidal  velocity  fields  ^  and  u  satisfy 
(2.3) 


(a) 

^(2o't-s)  -  ^Sf- 

2o<°> 

(b) 

=  -||  - 

The  linearized  perturbation  of  the  extra  stress  is  given  by 
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(2.4) 

(2.5) 


So  -  1  [So<»>]  >  So  -  »5oTt2o  -  i 

s=0 

OO 

S  =  Ji^Go(s)  I  3g]  ,  3G  =  V^V}^  +  V^TV£ 

S“0 


where  3^  is  the  first  functional  derivative  of  3",  evaluated  on 
the  basic  flow.  That  is,  3^  is  a  linear  operator  actiong  on  6G. 
We  assume  an  integral  representation  for  3^: 


09 

(2.6)  q  =  |  Kts^JdGtsJds 


The  isotropy  of  3F  implies  that 


(2.7)  2  So  ST  =  ?  [s  So  ST]],0 

(2.8)  2  5.  aT  -  -i.  [a  So  2Tl2  «S  2f]s=o 


which  holds  for  all  orthogonal  £  and  symmetric  <SG. 

For  simple  shear,  with  shear  rate  k  *  we  have 

n 


(2.9) 

Xq  ■  (I  -  sN)x 

(2.10) 

■  -s(N  +  NT)  +  s2NTJj[ 

(2.11) 

6G  =  V£  +  V£T  -  s(V£TN  + 

where 

(2.12) 

/o  1  o\ 

N  =  K  10  0  0 

,0  0  0 


The  extra  stress  of  the  basic  shear  flow  is  given  by  the  viscometric 
functions 


where  t  is  the  shear  stress  and  =  T11  -  T 22 ,  N2  =  T22  -  T33 
are  the  first  and  second  normal  stress  differences. 

Consistency  relations  between  the  kernel  function  K  and  t 
viscometric  functions  can  be  obtained  as  follows.  First,  we 
consider  families  g(A)  of  orthogonal  matrices  such  that  2(0)  = 
and  then  differentiate  (2.7)  with  respect  to  A,  at  A  =  0.  Sin<__ 
2(0)  is  an  arbitrary  skew-symmetric  matrix,  this  will  yield  three 
matrix  equations.  A  fourth  matrix  equation  is  obtained  by 
differentiating  (2.7)  with  respect  to  k,  i.e.  by  perturbing  the 
shear  rate.  Using  the  notation  (K) ^ j £m  =  k^, 
following  four  equations  [11] 


we  find  the 


U.15) 


-2t  N 


N, 


1 

2t 


(2.16) 


N!  +  N2 


0 

0 

0 


00 

=  (  2<sj^Kn  -  K22  -  KSK^J 


ds 


N1  +  N2 


=  |  -2icsK23ds 


(2.17) 


0  0 

0  0 

T  N. 


N. 


' 

=  2ks(ksK23  -  K13)ds 


Clearly,  only  (2.14)  and  (2.15)  are  relevant  for  two-dimensio..al 
flows. 

Equation  (2.8)  for  the  isotropy  of  -3^  can  be  used  to  show 
[11],  by  choosing  special  forms  for  Q,  that 


<2‘18)  Kij Jim  -  0 

whenever  there  are  an  odd  number  of  3's  in  the  set  {i,j,Jl,m}. 
Finally  the  symmetry  of  S  implies  that 

(2-19)  Ki j Jim  *  KjiJlm 


and  the  incompressibility  condition  implies  [19] 


K 


(2.20)  I  Kl!L  +  2ksK12  +  ksKi;l  =  0 


Prom  (2.6)  and  (2.11),  we  find  the  following  expression  for  the 
perturbed  extra  stress 


(2-21>  2=1  1,6=1  Sa6  SGc6  as 


where 


2  j*  =2^1,1  +  =1(^1,2  +  ?2,1}  +  =13^1,3  +  ^3,: 

*  A 


+=23(?2,3  +  ?3,2  "  ks?1,3)  +  (=33  "  =22^3, 3}ds 


(2.22) 


K,  =  -Ks  K22  +  K12 


(2.23) 


&2  “  =11  “  =22  ”  ks=12 


Note  that  K^,  k2  ,  and  K^,  are  all  related  to  viscometric 

functions  by  the  consistency  relations  (2.14)  -  (2.17),  but  that 
-  K,,  does  not  occur  in  any  of  the  consistency  relations. 


k 
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3 .  SHORT  MEMORY 


We  consider  in  this  section  the  short  memory  assumption  as 
applied  by  Akbay,  Becker,  Krozer  and  Sponagel.  In  general, 
smooth  deformation  histories  may  be  approximated  asymptotically 
in  terms  of  the  Rivlin-Ericksen  tensors. 


(3.1)  G( s)  -  1  ^[ult)]  +  0 


(SN+1) 


where 


A  fu(t)1  =  (-1) n  G(s) 


is  the  nth  Rivlin-Ericksen  tensor.  For  the  perturbation  term  6G, 
we  obtain 

(3.2)  6G  =  -s  A^uCt)]  +  0(s2) 

=  -2s  fi[u(t,]  +  0(s2) 

1  T 

where  D[u]  =  j(Vu  +  Vu  ) .  Under  the  assumption  that  this  leading 
term  is  sufficient  for  a  short  memory  fluid,  (2.21)  may  be  written 


(3.3) 


=  -2I,s1Li  ^  D°s}ds 

f°° 

-  “2 1  s^n  ~  S22)ds  35TT 

‘  o  i 

r°°  3u.  3u» 

-2  =Si2ds  ^  + 

'  o  z  i 

f  ^U3 

"2  sK13ds  <377  +  Jx7> 
o  ->  1 


*  +  t 
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3u-  3u. 

”2J  o  S“23ds  (3SJ  +  &q] 

00 

f  3u3 

-2  s(K33  -  K22)ds  ^ 

The  first  four  integrals  are  precisely  those  occurring  in  the 
consistency  relations/  and  if  terms  of  order  <s  can  be  neglected 
compared  to  others,  they  are  approximately  equal  to  expressions 
(2.14)  -  (2.17)  involving  viscometric  functions.  With  this 
approximation,  we  arrive  at  the  following  expression  for  the 
perturbed  extra  stress. 
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This  is  the  form  used  by  Akbay  et  al.  [ 2  J .  We  see  that  even  in 
this  approximation  there  remains  a  term  not  related  to  viscometric 
functions  (  see  [20]).  More  such  terms  would  occur,  if  higher 
order  terms  were  included. 

In  two-dimensional  flows,  (3.4)  reduces  to 


(3.5) 


-2t  Nx 
Nx  2t 


+ 


( 


3u.  3u_ 

_ £  +  _ —) 

ax2  +  sx^ 


The  equation  of  motion,  after  linearizing  at  plane  Couette  flow  and 
eliminiating  the  pressure,  becomes 


(3.6) 


p(3t  +  KX2  3x,}  ( 


3Uj 

3x. 


3u. 

4 

3x: 


12 


+ 


a2(Sll  -  S22> 
3x13x2 


We  introduce  a  stream  function  Hf  such  that  (u^,u2)  = 


Then  we  obtain 
(3.7) 


„  (9_  +  KX  2 _ i  +  iliL, 

P<3t  2  3xi  ^  ax/ 


=  (N  '  -  -i)  —  L¥  +  T '  IT4 '  +  "  — 


’1  k  '  3xl3m2 


K 


3 2  a2 

where  L  =  - *■  -  - 

3x2  3x1 


1  * 


,  3Y  3 jf  . 
l-3x2'”5x^}  * 


T 

2 
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4.  WHAT  IS  A  CHANGE  OF  TYPE  AND  WHY  SHOULD  WE  CARE  ABOUT  IT? 


We  shall  see  below  that  equation  (3.7)  will  have  an  instability 

if  a  certain  "Weissenberg  number"  exceeds  a  critical  value.  THis 

instability  will  be  associated  with  what  is  called  a  change  of 

type.  The  reason  why  this  is  important  is  that  change-of-type 

instabilities  have  quite  different  mathematical  properties  and  will 

lead  to  quite  different  physical  effects  whan  other  instabilities. 

In  order  to  explain  what  we  mean  by  a  change  of  type,  we 

must  first  explain  what  we  mean  by  the  symbol  of  a  differential 

operator.  To  obtain  the  symbol,  one  simply  replaces  the  derivatives 
a 

gj—  by  numbers  D^ .  Thus  we  can  assign  to  any  differential  operator 
a  polynomial.  For  example,  the  differential  operator  on  the  right 
of  (3.7)  would  have  the  symbol 

(4.1)  (N^  -  ^)d1d2(D22“D12)  +  t’<d22  “  di2)2  +  ]r  d12°2 2 

Two  differential  operators  are  said  to  be  of  the  same  type,  if  — 
up  to  a  transformation  of  the  independent  variables  —  their  symbols 
have  the  same  asymptotic  behavior  at  ».  If  the  asymptotic 
behavior  of  the  symbol  changes,  we  say  that  a  change  of  type 

a2  a2 

occurs.  For  example,  the  operator  — _  +  y— -  changes  type  as  y 

3x^  3y 

2  2 

goes  through  zero:  For  y  >  0  the  symbol  D1  +  yD2  is  a  positive 
definite  quadratic  form  (the  operator  is  then  called  "elliptic") , 
for  y  <  0  the  symbol  is  an  indefinite  quadratic  form  (the  operator 
is  then  called  "hyperbolic") . 

Mathematical  problems  such  as  existence,  uniqueness  and 
regularity  of  solutions  or  continuous  dependence  on  the  data 
depend  essentially  on  the  type  of  the  differential  operator.  For 
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second  order  operators  it  is  well  known  that  in  the  elliptic  case 
the  "good"  problem  to  look  at  is  a  boundary  value  problem,  while 
in  the  hyperbolic  case  it  is  an  initial  value  problem. 

Heuristically,  the  reason  for  this  is  that  the  type  describes 
how  the  operator  acts  on  rapidly  varying  functions.  For  such 
cases,  one  expects  that  locally  the  variation  of  the  coefficients 
can  be  neglected,  and  that,  in  the  interior  of  the  domain  where  the 
problem  is  posed,  the  influence  of  boundary  solutions  can  also 
be  neglected.  For  constant  coefficient  problems  without  boundaries, 
the  symbol  is  just  the  Fourier  transform  of  differential  operator. 


How  does  the  foregoing  apply  to  equation  (3.7)?  If  we  formally 


set  ft  * 


=  D1  =  ia. 


=  d2  =  i&. 


then  the  left  side  of  the  equation  becomes 

(4 .2)  p(0r+  Kx2ia)  (— B2-a2)  , 


and  the  right  side  becomes 

(4.3)  -(N^  -  ^)  (a2  -B2)a6  +  f(a2  -  B2)  2  +  ~  a2B2 


The  expression  (4.3)  is  a  homogeneous  fourth  degree  polymonial. 

It  turns  out  that  this  polynomial  is  positive  definite  if  and  only  if 

N  2 

(4.4)  We2  =  K2/(n(K)T'(K))  <  16 

2 

If  We  >  16,  the  polynomial  is  indefinite-  If  we  now  formally  put 

(4.2)  equal  to  (4.3),  then  we  find  that  a  is  negative  if  a,B  are 
2  2 

large  and  We  <  16.  If  We  >  16,  then  there  is  a  sector  in  the 
a, &-  plane  for  which  a  becomes  positive  for  a,B  large,  moreover 
o  becomes  arbitrarily  large  in  this  case  as  a,B  -*■  00 . 
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o 

This  heuristic  argument  suggests  the  following:  If  We  <  16, 

then  at  least  highly  oscillatory  modes  will  be  stable.  If 

instabilities  exist,  then  they  will  only  affect  a  finite  number  of 

2 

modes.  If,  on  the  other  hand.  We  >16,  then  we  expect  to  find 

an  infinite  number  of  unstable  modes.  In  addition,  the  growth 

constants  of  these  modes  can  be  arbitrarily  large.  Thus  the 

change  of  type  leads  to  an  instability  which  is  in  a  sense  much 

stronger  than  those  studied  in  bifurcation  theory.  One 

cannot  expect  a  secondary  flow  whos^  dynamics  would  be  governed  by 

the  evolution  of  one  or  several  modes.  Rather,  the  initial  value 

2 

problem  becomes  ill-posed,  i.e.  there  are  flow  fields  for  We  >16, 
which  would  not  occur  even  as  transient  states,  since  "random" 
disturbances  containing  all  modes  would  blow  up  instantly.  In 
situations  where  this  "forbidden"  unstable  region  lies  between  two 
stable  ones,  a  phase  transition  can  occur.  That  is,  there  can  be 
discontinuous  solutions  whose  values  lie  in  both  of  the  stable 
regions.  The  two  regions  are  then  referred  to  as  different 
"phases".  If  the  forbidden  region  stretches  all  the  way  to 
infinity  and  a  second  phase  does  not  exist,  the  result  of  a  change 
of  type  may  be  something  like  fracture. 
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5.  INSTABILITY  OF  PLANE  COUETTE  FLOW  WITH  SHORT  MEMORY . 


In  this  chapter,  we  give  a  more  detailed  analysis  of  equation 
(3.7),  and  we  provide  a  rigorous  basis  for  some  of  the  heuristic 
conclusions  obtained  in  the  last  chapter.  The  stream  function  Y 
is  represented  by  normal  modes 

(5.1)  Y  -  eot  eiaxl<|>(x2)  . 


We  wish  to  find  eigenvalues  a  and  eigenfunctions  $(*)•  For  simplicity, 
we  set 


(5.2) 


E  -  »  r  = 


FTkT 


( 


Ni(k)  • 


R  -  xtfk'  0  =  “iaK:  c- 


Moreover,  we  assume  that  the  width  h  of  the  channel  is  1. 

Equations  (3.7),  (3.8)  then  become 

(5.3)  iaR(x2-c)  ($'  '  -  ot24»)  =  4>1V  +  a2(2-E)<J>''  +  a4«|)  +  iar(4>"'  +  a2*') 

(5.4)  <P  =  <P'  ~  0  at  x2  =  0,  x2  =  1 

Let  us  assume  that  T'  >  0  euid  hence  R  >  0.  In  this  case  we  want 
to  show  that  when  T2  <  4E  the  imaginary  part  of  c  must  be  negative 
provided  that  a  is  large  enough.  For  this,  we  form  an  energy 
integral 

i 

f  [$iv  +  a2 (2-E)  <j>*  '  +  a4<J>  +  iar(<j>'"  +  a2*))'  )  ] ^dx2 
*  o 

i 

=  iaR  [  (x2-c)  ( 4> *  '  -  a2<j>)<Fdx2 


(5.5) 
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and  integrate  by  parts#  to  find 


[  |tj>"  +  a2<j>|2  dx2  +  Ea2  [  |<j>'|2  ^2  “  iar  f  +  a2<p)<p'  dx 

■*0  ^  0  '  0 

1  1 

(5.6)  =  -iaR  I  (x2-c)  [|<j>'  |2  +  a2|4>j2]  dx2  -  iaR  f  <J>'$dx2. 

•  A  *  A 


There  is  a  constant  y  such  that 

o 


Hr 


*a*2  i  Yo 


f  If '  I 

'  n 


’dx2  for  any 


♦  satisfying  $(0)  =  <|>(1)  =0.  Taking  the  real  part  of  (5.6)  and  using 

this  inequality,  we  find 

1 

(5.7)  -aR  Imc  [  (|<J>’|2  +  ot2 1 4>  |  2)  dx2 

1 

>|  |  <j) '  '  +  a2<|>  |  2dx2  +  (E - ^-)  a2  j  |  <j> '  |  2dx2 


+  af  In  f  ( *  *  +  a2$)tf>'dx. 


By  applying  Cauchy-Schwarz  to  the  last  term,  we  find  that  the 

2  RYo 

right  hand  side  of  this  is  positive  as  long  as  r  <  4(E  -  ) .  Thus, 

2 

if  r  <  4E  and  a  is  large  enough,  we  must  have  Imc  <0.  If  aRx2 

is  put  to  zero  in  (5.3),  the  term  proportional  to  R  in  (5.7) 

2 

disappears,  and  stability  for  any  01  is  guaranteed  if  r  <  4E. 

In  the  following,  we  neglect  the  Rx2  term.  Then  with 
X  *  -iaRc,  we  have  (5.3)  in  the  form 

(5.8)  AU"  -  a2<p)  =  <J>1V  +  ( 2-E)  a2<f> '  '  +  a4<{>  +  iar(<f>"'  +  a2<J>'). 


We  shall  prove  that  there  is  an  eigenvalue  X  *  0  at  a 
value  of  IMa)  which  converges  to  /3E  as  a  00  (This  proof  is 
due  to  M.  Ahrens).  Since  (5.8  )  is  a  linear  fourth  order  ODE, 
it  can  be  solved  exactly.  We  look  for  exponential  solutions  of 
the  form  <p  (x2)  =  elarx2  and  obtain  the  following  algebraic  equation  for 
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(5.9)  r4  +  Tr3  -  (2-E)r2  -  Tr  +  1  +  ^-(r2  +  1)  =  0 

ex 

We  are  interested  in  neutral  stability, X  =  0,  so  that 

(5.10)  r4  +  Tr3  -  ( 2-E) r2  -  Tr  +  1  =  0 

Equation  (5.10)  is  a  fourth  degree  polynomial  in  r  and  hence 

r2 

solvable  by  standard  methods.  Letting  Y  =  4 -  E,  80  that  y  <  0 

corresponds  to  linear  stability;  «/e  can  write  the  discriminant 
for  (5.10)  in  the  form 

(5.11)  A  =  16y2{16y  +  (E+4)2}. 

Thus  A  >  0  for  y  >  0,  assuring  us  that  (5.10)  has  distinct 
roots  in  the  case  of  interest.  In  addition,  the  four  roots  of 

(5.10)  are  real.  For  completeness,  we  list  the  roots  here. 

,  k  =  1,  2,  3,  4  where  rfc  are  the  roots  of 

(5.10)  .  Thus 

(5.12)  Rx  =  +  ^T2  +  /=yj} 


r 

Define  j 


where  y^,  y2,  Y 3  are  real  and  negative  and  are  solutions  of  the 
cubic 

(5.13)  Y3  +  bjY2  +  b2Y  +  b3  =  0, 
where 
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where 

bx  »  2(|r2  +  2-E) 

(S' *1^)  b2  *  IS  P4  +  (i-E)1,2  “  4E  +  E2 

b3 =  y2  v2-™2 

It  is  easy  to  check  that  b1#  b2,  b3  are  all  positive  for  the  case 
of  interest,  r2  >  4E,  if  we  assume  E  >  -4.  The  rule  of  signs 
applied  to  (5.13)  thus  shows  y1#  y2,  y3  are  negative,  so  that  rk, 
k  =  1,  2,  3,  4  are  all  real,  as  asserted  previously.  In  fact,  the 

y^,  ^2 '  ^3  ma^  be  f°un<*  f rora 

(5.15)  y1  m  +  YI  +  ITT} 

y2  -  |-{-b1  +(-|  -  ^  un  *  (-1  +  ^ 

Yj  «  y{— +  ^  -y  i)  ¥l  +  l"y  —  i)  ^XI } 

where 

(5.16)  I  =  IT  =  -(E  +  4)  (18y  +  (E  +  4)2)  +  |/I/A  i. 


We  are  now  able  to  write 
type  problem,  in  the  form 


(5.17) 


f 


4 

Z 

k=l 


ak  e 


iotrkx2 


down  a  solution  to  the  Orr-Sommerfeld 


where  r^  are  determined  by  (5.12),  (5.13)  above,  and  the  a^  are 
constants.  The  boundary  conditions  allow  a  non-trivial  <j>  (x2)  if  and 
only  if 
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This  implicitly  defines  a  relationship  T(a)  which  must  hold  for 
neutral  stability.  Reduction  to  the  yk,  k=l,  2,  3  given  in  (5.15) 
gives  (5.18)  in  the  form 

(5.19)  f(a,f)  =  (y3  -  y2)  cos  (oti^y^)  +  (yx  -  y3)  cos  (rf/^yj) 

+  (y2  -  Y1)cos(o</^J)  =  0 


Recall  that  y^,  y2#  y3  are  independent  of  a,  and  are  smooth 

functions  of  r  in  our  case  of  interest,  r  >  4e. 

Equation  (5.19)  was  numerically  solved  to  display  the 

neutral  stability  curve  T(a),  for  several  choices  of  the  parameter 

E.  A  typical  case,  E  =  2,  is  given  in  Figure  1.  Note  the  behavior 

2 

of  T  as  a  +  *,  indicating  that  T  =  4E  is  a  horizontal  asymptote 

2 

We  shall  show  that  P  =  4E  is  indeed  a  horizontal  asymptote. 

2 

Thus,  loss  of  linear  stability  occurs  at  this  value  of  r  for 
disturbances  with  short  wavelength  (a  large) .  The  numerical 
results  of  Akbay  and  Sponagel  [4]  indicate  the  same  result. 
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Asymptotic  Behavior 


We  will  take  the  limit  of  (5.19)  as  T  4-  4E,  i.e.  as  y  4-  0, 
First,  from  the  form  of  y^  y 2,  y^  given  in  (5.15)  we  find 


(5.20)  y1  +  -(E  +  4) 


y2'  y3  -  0 


as  y  +  0. 


a  jl 

Then  computing  f  f—  from  (5.18),  and  using  this  last  result  (5.20) 

<7  Qt  / 


gives 


(5.21)  f,  |f  -  0  as  y-  0; 


We  now  consider  —  .  From  the  definition  of  y^,  y2,  y^,  we  can 


compute 


1  _  E  +  8 


(5*22) 


E  +  4 


as  y  4-  0. 


Hence,  (5.19)  implies 


( 5 ’ 23 )  -E  /E{cos (a/E  +  4)  -  1  +  -  -2-—  a2}  as  y  +  0. 


g(a)  =  {cos(a/E  +  4)  -  1  +  —  ^  4  a2}. 


Then 


g(0)  =  0  and  g 1  ( ot)  =  Ze  +  J{ZE  +  4a  -  sin(a/E  +  4)  } . 


Clearly  g'(a)  >  0  for  a  >  0.  Thus,  g(a)  >  0  for  a  >  0. 

g  f 

Hence,  for  a  >  0,  does  not  reduce  to  zero  as  y  +  0. 

The  behavior  of  f  near  y  =  0  is  approximately  described  by 


(5.24)  f  (a,  T)  *  y^costaZ-yj)  -  cos  (a/-y3> ) 

3y2  3y3 

Since  yy-  and  yy-  have  different  limits  as  y  +  0, 

the  term  in  brackets  will  have  an  increasing  number  of  zeros  for 

small  values  of  y  as  a  grows  large.  It  is  clear  that  this  property 

persists  for  the  exact  f.  Hence  we  find  an  infinite  number  of 

2 

curves  r  =  T(a) ,  which  must  approach  the  line  r  =  4E  as  a  +  ®. 
Since 


(5.25) 


dr  _  ,3f  -1  3f 

3a  “l3r'  3a' 


the  slope  of  each  of  these  curves  must  approach  0  as  a  “.  Hence 
2 

T  =  4E  is  a  horizontal  asymptote  for  these  curves. 

We  turn  next  to  a  description  of  results  obtained  by  numerical 
calculation  at  finite  Reynolds  numbers.  Figure  2  shows  the 
neutral  stability  curve  for  the  case  R=l,  E-4,  showing  an 
asymptotic  behavior  r  *►  4  as  a  <».  The  method  of  calculation 
follows  that  of  Orszag  [17] .  Equation  (5.3)  is  discretised  by 
representing  <j>  as  a  sum  of  Chebyscher  polynomials.  The  resulting 
matrix  equation  is  solved  for  c,  given  a  and  T  .  The  computations 
were  checked  against  Gallagher  and  Mercer's  [12]  results  for  the 
case  r  =  0,  R  =  1,  10,  20,  30  using  20  Chebyshev  polynomials. 

The  number  of  chebyshev  polynomials  used  for  Figure  2  varies  from  40 
at  o  =  20  up  to  100  at  a  =  50.  For  negative  r,  the  neutral  stability 
curve  is  a  reflection  across  the  a-axis.  When  R  is  increased 


Figure 


this  symmetry  is  lost  and  the  neutral  stability  curve  for  r  positive 
is  raised.  These  trends  are  shown  for  a  -  4  in  Figure  3 
for  which  40  Chebyshev  polynomials  are  sufficient  for  4-digit 
accuracy  in  c.  A  similar  calculation  at  a  =  50  requires  over 
80  Chebyshev  polynomials  for  the  same  accuracy.  At  R  =  1 
the  neutral  stability  curve  lies  near  r  =  4.00  while  at 
R  *  25,  it  lies  near  T  =  4.04.  Hence  it  appears  that  the  curve 
asymptotes  to  r  =  4  as  a  ■*  °°  but  more  slowly  as  R  increases. 
Curiously,  the  influence  of  finite  R  seems  to  be  stabilizing 
rather  than  destabilizing,  at  least  for  low  values  of  R. 
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FIGURE  3 


TABLE  OF  THE  WEISSENBERG  NUMBERS  FOR  WHICH  Re  O  *  0 
a  ■  4  FOR  VARIOUS  REYNOLDS  NUMBERS 


R 

jf 

Positive  — 
n 

N,  , 

<r> 

V 

Negative  — 

Ni  • 

(r) 

1 

4.4 

-4.4 

5 

4.5 

-4.4 

15 

4.7 

-4.4 

K 
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6. _ A  CRITICAL  ASSESSMENT  OF  THE  SHORT  MEMORY  THEORY 


A  difficulty  in  the  theory  of  Akbay  et  al.  is  that,  if  the 
short  memory  approximation  is  applied  to  the  basic  flow  one  finds 
that  the  first  normal  stress  is  negligible  and  hence  P  »0.  That 
is,  if  short  memory  is  assumed,  then  all  terms  in  (3.1)  except 
the  first  can  be  neglected,  and  we  get  the  Reiner-Rivlin  fluid, 
which  has  N1  b  0.  The  instability  criterion  will  not  be 
achieved.  We  illustrate  this  for  the  following  model  constitutive 
law  of  Oldroyd  type 

(6.1)  £  +  A  IS  +  (v*7)S  +  U’S  -  S  *oj)  +  a(D«S  +  S-D)  ]  =2yD 

where  S  is  the  stress  and  0  and  are  the  symmetric  and  antisymmetric 

parts  of  the  velocity  gradients,  y  and  X  are  arbitrary  positive 
constant,  and  a  is  supposed to  be  between  +1  and  -1. 

In  simple  shear  flow,  we  get 

T  * _ yjs _ 

1  +  k2X2 (l-a2) 


N1  * 


2yAtc _ 

2  2  5 

1  +  iczAz(l-az) 


If  |a|  <  1,  t'  becomes  negative  for  large  k.  This  will  lead  to  an 

instability,  but  not  of  the  kind  considered  by  Akbay  et  al. 

2  2 

With  n  =  (1  -  a  ) X  ,  it  can  be  shown  that 


16  (  K  } 


X2k2(1  - 


<2n) 


,  2  2 
X  K 


TT 


4  (1  +  k  n) 


Hence  the  instability  criterion  does  not  apply  if  X<  is  small. 
The  maximum  value  of  the  left  hand  side  is 
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- i - 5-  •  (/Z  -  l)2, 

4(1  -  a^) 

which  never  reaches  1  unless  a  is  sufficiently  close  to  1. 

The  problem  of  consistency  in  the  theory  of  Akbay  et  al.  raises 

the  question  whether  changes  of  type  (or  their  absence)  can  be 

discussed  without  resorting  to  a  short  memory  hypothesis. 

For  a  K-B  K  Z  fluid  [15]  ,  [8],  this  has  in  fact  been  done  [22]. 

The  K-B  K  Z  model  includes  the  special  cases  a  =  ±  1  of  (6.1). 

In  these  cases*  the  condition  (1.3)  is  satisfied  for  large  enough 

K,  however,  it  was  shown  in  [22]  that  a  change  of  type  does  not 

occur  I  Thus  the  criterion  (1.3)  is  invalid  without  the  short 

memory  hypothesis.  Thus  the  question  whether  large  values  of 

would  produce  instabilities  in  shear  flow  remains  open. 
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